Critical dynamics of a uniaxial and dipolar ferromagnet 
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We study the critical dynamics of three-dimensional ferromagnets with uniaxial anisotropy by 
taking into account exchange and dipole-dipole interaction. The dynamic spin correlation functions 
and the transport coefficients are calculated within a mode coupling theory. It is found that the 
crossover scenario is determined by the subtle interplay between three length scales: the correlation 
length, the dipolar and uniaxial wave vector. We compare our theoretical findings with hyperfine 
interaction experiments on Gd and find quantitative agreement. This analysis allows us to identify 
the universality class for Gd. It also turns out that the /iSR relaxation rate can be best fitted if it 
is assumed that muons occupy octahedral interstitials sites within the Gd lattice. 
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I. INTRODUCTION 

The spin dynamics of simple ferromagnets in the vicin- 
ity to their Curie point is an archetypical example for 
critical dynamic phenomena near second-order phase 
transitions. In recent years it became increasingly clear 
that the dipolar interaction has a dramatic efFect_on the 
critical spin fluctuations of all real ferromagnets 
mode coupling theory (for a review see e.g. Ref. 
eluding the dipolar coupling has lead to remarkable suc- 
cess in explaining the experimental data for materials 
with cubic lattice structure (e.g. Fe, Ni, EuO, and EuS). 

In this paper we extend this analysis to anisotropic 
magnetic systems. A short account of part of these re- 
sults with emphasis on the interpretation of experiments 
on Gd has been given recently [^. Anisotropy acts to 
suppress critical fluctuations perpendicular to the easy 
axis of magnetization and breaks certain local conser- 
vation laws. This has a marked effect on the critical 
dynamics of spin fluctuations near Tc- Combining the 
effect resulting from magneto-crystalline anisotropy and 
long-ranged dipolar interaction one expects a crossover 
scaling behavior with three scaling variables. Besides the 
correlation length ^ there are two more length scales re- 
sulting from the strength of the dipolar interaction and 
the anisotropy energy with respect to the exchange en- 
ergy 

Dipolar interaction is known to be a relevant pertur- 
bation with respect to the fixed point of a rt-component 
Heisenberg model, where the spins are coupled by short- 
range exchange interaction. It drives the system to a new 
dipolar fixed point, whose nature dramatically depends 
on the number n of components of the order parameter. 
For isotropic Heisenberg ferromagnets (n = 3) the re- 
sulting isotropic dipolar fixed point is characterized by a 
set of critical exponents which are only slightly different 
from the corresponding values at the (isotropic) Heisen- 
berg fixed point. The consequences of this crossover on 



the static and dynamic correlation functions is by now 
well known and has been reviewed recently For uni- 
axial {n — 1) systems, however, it was shown by Larkin 
and Khmelnitskii j?) that dipolar interaction asymptoti- 
cally leads to classical critical behavior with logarithmic 
corrections in three dimensions. The asymptotic behav- 
ior of this system was also studied by means of renor- 
malization group theory ||^-[lO|], which revealed that the 
one loop calculation agrees with the asymptotic results 
of Ref. 0. The complete crossover from Ising behav- 
ior with non classical exponents to asymptotic uniaxial 
dipolar behavior has also recently been analyzed within 
a generalized minimal subtraction method . 

In anisotropic ferromagnets with an easy axis 
anisotropy there are two relevant perturbations, the dipo- 
lar coupling g and the anisotropy parameter m. As a 
consequence upon approaching the critical temperature 
the system passes through a considerably more complex 
crossover region before it reaches its asymptotic critical 
behavior. Four nontrivial fixed points determine the flow 
of the various model parameters: the isotropic Heisen- 
berg, the uniaxial Ising, the isotropic dipolar and the 
uniaxial dipolar fixed point. Depending on the relative 
strength of the dipolar and uniaxial anisotropy different 
scenarios are possible. For g > m the system is supposed 
to show a crossover cascade from isotropic Heisenberg to 
isotropic dipolar to uniaxial dipolar critical behavior. For 
m > g the system first crosses over from the Heisenberg 
to the Ising fixed point before it turns to the asymptoti- 
cally stable uniaxial dipolar fixed point. Both cases seem 
to be realized in nature; e.g. one finds m ~ 1.2 10~^ and 
g « 2.0 10-3 for FeuNdsB (see Ref. @, whereas m< g 
in Gd (see below)). 

In this paper we study the dynamics of such anisotropic 
ferromagnets under the combined infiuence of dipolar in- 
teraction and magnetocrystalline anisotropy. We proceed 
as follows. In Section || we discuss the model Hamilto- 
nian and the effect of dipolar interaction on the critical 
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behavior in cubic and hexagonal lattices. We define the 
dimensionless parameters characterizing the strength of 
the dipolar and the uniaxial anisotropy. In Sec. Ill 



we 



discuss the static critical behavior of uniaxial dipolar fer- 
romagnets and derive the eigenvalues and eigenvectors of 
the static susceptibility matrix. The mode coupling the- 



ory in terms of these eigenmodes is formulated in Sec. IV 



discussion of the results. Some of the technical details of 
the theory are collected in the appendix. 



II. THE MODEL OF AN UNIAXIAL DIPOLAR 
FERROMAGNET 



na/3 (gL^^s f J 47r f 3qaqp 



(2.4) 



We briefly discuss the limits of an isotropic Heisenberg 
magnet without dipolar interaction and of an isotropic 
Heisenberg magnet including dipolar interaction. The 
general case of an uniaxial and dipolar ferromagnet is 
discussed in detail. The mode coupling equations are 
solved analytically in certain limiting regions of param- 
eter space and analytically for intermediate parameter 
values. The resulting crossover scenarios for dominating 
magneto-crystalline and dipolar anisotropy are discussed, 
respectively. In Sec. we compare our theoretical find- 
ings with results from vari ous hyperfine interaction ex- 
periments on Gd. In Sec. VI we give a summary and ^ 



where Va is the volume of the primitive unit cell with lat- 
tice constant a, and are constants, which depend on 
the lattice structure (see e.g. ||^). Upon expanding the 
exchange interaction. 



Jq = Y^.h^e'^'^' -Jo- Jq^a" + 0{q^) , (2.5) 

i 

and keeping only those terms, which are relevant in the 
sense of renormalization-group theory this results in the 
following effective Hamiltonian for dipolar ferromagnets 
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q 



(-Jo + Jq^a^) 6ai3 + Jg 



gag 13 

g2 



5" 5^, (2.6) 



where the Fourier-transform of the spin variables is de- 
fined by 



(2.7) 



We consider a system with N identical spins fixed on 
the sites of a three dimensional lattice. Taking into ac- 
count magneto-crystalline anisotropy as well as dipolar 
interaction it is described by a Heisenberg hamiltonian 



H 



Blasts'; 



(2.1) 



The magnitude of the magneto-crystalline anisotropy of 
the system is given by A = jH/J^. Here we focus on 
uniaxial anisotropy, A > 1, with the easy axis of mag- 
netization along the z-axis. The dipolar interaction is 
characterized by the tensor 



Here we have defined a dimensionless quantity g as the 
ratio of the dipolar energy (gLMB)^/^'^ and exchange en- 
ergy 2 J, multiplied by a factor Ana'^/va, which depends 
on the lattice structure. 

Aira'^ [gL^iBY I a? Dipolar Energy 
g = — oc 



2J 



Exchange Energy 



Strictly speaking there are dipolar corrections of order 
0{q^) to the exchange coupling. But, those can be ne- 
glected, since the strength of the dipolar interaction is 
small compared with the exchange interaction. 

For Bravais lattices with a hexagonal- closed packed 
(hep) structure the dipolar tensor to leading order in q 
becomes of the form (see Rcf. llB] and the appendix) 



_ {gLt^Bf 
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Qaqp 
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(2.2) 



with Xij ~ Xj — Xj, gL the Lande factor, and /is the Bohr 
magneton. As shown by Cohen and Keffer jl3| dipolar 
lattice sums 



(2.3) 



can be evaluated by using Ewald's method |lj]. For infi- 
nite three-dimensional cubic lattices one finds to leading 
order in the wave vector q |15[ 



(2.9) 



with the coefficients given in table ||. In the same way as 
above we find 



- Jo" - Jm" + Jg^a^ ]Sap 



J qaqp 



qa qf3 

-q q ■ 
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where the contribution of the dipolar interaction to the 
uniaxial anisotropy is given by 
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(2.11) 



Note that the term proportional to qaQp/q'^ depends on 
the lattice structure only via the volume o f the unit cell. 
Therefore the value of g is identical to Eq. |2.8| . 

There are two s ource s of uniaxial anisotropy in the 
Hamiltonian, Eq. 2.10, magnetocrystalline anisotropy 
Jg", defined by J° « - Jq^a"^ + C (<?"*), and m" which 
characterizes the dipole-dipole interaction. Crystal field 
contributions are expected to be the dominant factor in 
systems like LiTbF4 ||l|l and Fe^NdaB ||. 

In addition, the dipolar interaction introduces an 
anisotropy of the spin-fluctuations with respect to the 
wave vector q which is reflected by the term propor- 
tional to qaqjs/q^- The magnitude g of this anisotropy 
is given hy g = 47r(gi/is)^/2 Jwa- We define a dimen- 
sionless quantity 



f3i)/2Jva 



(2.12) 



proportional to the ratio between the anisotropy energy 
and the exchange energy. Putting in values for Gd the 
ratio of the dipolar contribution to the term qaqp/q^ and 
to the uniaxial anisotropy is yjg/m = 7.8738. In sec- 
tion we will show that all available data for Gd can 
be explained by assuming that the uniaxial anisotropy is 
solely due to the dipolar interaction. 



III. THE CRITICAL STATIC BEHAVIOR 



A. The static susceptibility 

Upon using standard techniques such as the Hubbard- 
Stratonovich transformation one may derive an effective 
Landau-Ginzburg free energy functional from the micro- 
scopic Hamiltonian, Eq. ( p.l[ ). Then the renormalized 
free energy in Gaussian approximation becomes, 

^=^x;^(k)5"(k)5'^(-k), (3.1) 

where = Va J d^k/(2Tr)^. The inverse renormalized 
propagator (susceptibility) is given by 



[ra + a k )6ai3+g 



(3.2) 



defines the correlation length i = £.o[{T -T^) /Tc]'" with 
a non-universal amplitude and the effective exponent 
V. In the hard sector the "masses" do not vanish upon 
approaching the critical temperature, but saturate at a 
finite value m 



(3.4) 



The parameter m = (g^a)^ characterizes the magne- 
tocrystalline anisotropy. The ratio of dipolar to ex- 
change interaction can be described by the Parameter 
g — [q^a)"^. Thus, in addition to the correlation length 
there are two more relevant length scales, q~^ and q~^- 
With the scahng variables, R = (x, y, z), 



(3.5) 



1 (Id 

qi q 

the susceptibility matrix becomes, 



Xa/3(k) = Xa/3(k;^,go,9^) = J (ag) ^Xa/3(k;R), (3.6) 



with 



Xafj 




y''-^ (3.7) 



For a more complete discussion of the static crossover 
and a calculation of the effective exponents we refer the 
reader to the literature [n]jl|-l2^. 



B. Eigenvalues and eigenvectors 



To investigate the physical properties of the static sus- 
ceptibility one has to find the eigen vect or s an d eigenval- 
ues of the spin system. ^From Eq. (3.2)-( p!4[ ) the eigen- 
values of the inverse susceptibility matrix are found to 
be 



where we have taken an Ornstein-Zernike functional form 
and assumed that the static crossover can be described in 
terms of scale dependent parameteres r^ and g. In other 
words, all the crossover is contained in effective static ex- 
ponents. We have chosen the reference frame such that 
the z-axis coincides with the easy axis of magnetization. 
Hence the "mass" of the corresponding spin fluctuations, 



C2 ' 



(3.3) 



Ai(k) = + 
X2{k)^k^ + r^ 



1a 

1 r 2 



w] 



X,{k)^k' + C' + -[ql+ql-W] 



where W — \J [q^, + q\Y - ^q^^q^k^ / k'^ . 
spending eigenvectors are given by 



(3.8) 
(3.9) 

(3.10) 
The corre- 
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ei(k) 



k2 

ki 



( 



62 (k) = 







y ^fcf+fci+/|(k) / 
fci 



sgn(fci) , 
\ 



(3.11) 



sgn(/+(k)), (3.12) 



63 (k) 



fc2 



/-(k) 



sgn(/_(k)), (3.13) 



with 



/±(k) 



2fc3(l-fc|/fc^)g^ 
g2 +(l_2fc|/fc2)q2 



(3.14) 



This eigenvectors obey the symmetry relation 6Q(k) = 
eQ(— k). It is interesting to note that due to the com- 
bined effect of the dipolar interaction and the uniaxial 
anisotropy the eigenvalues of the susceptibility matrix 
remain finite in the limit fc ^ and upon approaching 
the critical temperature. Only if the angle v between 
the easy axis of magnetization and the wave vector is 
V = 90° the third eigenvalue becomes critical. In the 
latter case the eigenvectors reduce to ei(k) = (^2, ki,0), 
e2(k) = {ki,k2,0) and 63(k) — (0,0,1), i.e. the critical 
eigenvector is along the easy axis of magnetization. 




FIG. 1. Eigenvectors in the uniaxial limit. The easy axis 
of magnetization points along the z-axis. The eigendirections 
in the uniaxial case are the easy axis of magnetization and 
two directions perpendicular to the easy axis, which here are 
chosen such that one is the projection of the wave vector into 
the xy-plane. 

In order to understand the critical behavior of the 
static susceptibilities it is useful to consider the uniaxial 
and the dipolar limit. If the dipolar anisotropy can be ne- 
glected with respect to the uniaxial anisotropy, g^-, <C , 
the eigenvectors reduce to (see figure |l|) 




sgn(fci 



sgn(fci 



61 (k) = 



e2(k) 



e3(k) = 



The corresponding eigenvalues are given by 



(3.15) 



(3.16) 



(3.17) 



Ai(fc) = fc2+^-2 
X3{k)^k^+r^. 



(3.18) 
(3.19) 
(3.20) 



In the pure dipolar limit, = 0, one finds the eigenvec- 
tors 



k2 



61 (k) 
e2(k) 

e3(k) 



\ 



\ / 

^sgn(fc3), 
/ fci ks 



sgn(A;i 



k2 k-j 



(3.21) 
(3.22) 

(3.23) 



where 62 (k) points along the wave vector k and 61 (k) is 
perpendicaular to the wave vector k and the easy axis of 
magnetization. The third eigenvector is perpendicular to 
k and ei. 




FIG. 2. Eigenvectors in the dipolar limit. The easy axis of 
magnetization points along the z-ajcis. The eigendirections in 
the dipolar case are paralell and perpendicular to the wave 
vector k. 
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The corresponding eigenvalues 



Ai(k) = fc2 
A2(k) - 
A3(k) = 



(transverse Ti), (3.24) 
(longitudinal L), (3.25) 
(transverse T2). (3.26) 



show that the spin fluctuations transverse to the wave 
vector k are critical, but the longitudinal ones stay finite. 
The dipolar wave vector determines the crossover scale 
where the longitudinal susceptibility turns finite. 



into account, the spin relaxation functions enter quadrat- 
ically into the coupled integro-differential equations for 
the Taj3 (q, t) . Frequently one introduces in addition a 
Lorentzian approximation for the Kubo relaxation func- 
tions, which results in a simplified set of mode coupling 
equations for the line widths. For instance the Resibois- 
Piette scaling function for isotropic fcrromagnets is ob- 
tained on this level of approximation |2l| . 



A. Equations of motion for the eigenmodes 



IV. MODE COUPLING THEORY FOR THE 
CRITICAL DYNAMICS 

In this section we calculate the wave vector and fre- 
quency dependence of the spin correlation functions, 
quantities which are (in principle) directly accessible to 
neutron scattering experiments. For this purpose we use 
standard mode coupling theory. The basic idea underly- 
ing MC theory is that near the critical point the relevant 
dynamics are described through slowly varying macro- 
scopic modes; i.e.: the conserved quantities and the or- 
der parameter. The dynamics are formulated most con- 
veniently in terms of the Kubo relaxation functions, 

00 

<i>„0(q,i) = zlmr j dTe-^^([s"(q, r), s^(q, 0)^]), (4.1) 
t 

for the components s" (q, t) of the eigenvectors of the 
static susceptibility matrix. Here we use the normaliza- 
tion <i>Q^(q, i = 0) = 1, i.e., the spin variables are nor- 
malized with respect to the static susceptibilities. Here 
( ) denotes the thermal average and [ , ] the commuta- 
tor. The corresponding frequency dependent relaxation 
functions are defined by a half sided Fourier transform 

00 

^o.p{ci,uj)^ j dte'^'<^^p{ci,t). (4.2) 


This Kubo relaxation function is related to the transport 
coefficients Tap{^., t) and the frequency matrix tt'Q,^(q) by 

^ $a/3(q, t) = iwa^(q)$^0(q, t) 

- f dTTa.^{q,t-T)^^f3iq,T), (4.3) 
Jo 

where the frequency matrix is given by 
1 



Wa/3(q) 



VXa(q)x,3(q) 



([s"(q),s'^(-q)]). (4.4) 



The coefficients, To.(3{ci,t), of the memory matrix are 
determined self-consistently from decay processes of the 
spin-modes. If only two mode decay processes are taken 



The analysis of the static susceptibility suggest to de- 
compose the spin operator in three components along the 
eigenvectors of the spin Hamiltonian, 



Sk = Skei(k) + 4e2(k) + sle^ik). 



(4.5) 



Written out in its cartesian components (which from now 
on are indicated by Latin labels) this reads 



5^ — s^eaiCk). 

a 

The back transform reads 



(4.6) 



(4.7) 



where the components of the spin operators in the eigen- 
vector basis are labeled by Greek indices. With the com- 
mutation relations for the cartesian components of the 
spin operator one is led to the commutators of the spin 
components s". 



1 

q^k 



(4.8) 



with 



f^a/sC*' = X! ^^J*^ e^i(q) e7fc(q - k), (4.9) 

ijk 

where Sijk is the Levi-Cevita symbol. The Heisenberg 
equations of motion become 



dt 



a 



} 



(4.10) 



In the classical limit this can be rewritten as 
= 2 J 

with 



^=2j/5:Tf''(k,q).(^.^_„ (4.11 
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- -^i(kjL/iQ(K, q) 






— — A2(q- 




Ta (k,q) 


= (-^i(k) - A3(q- 


^^^^ tt^ /i^ ^\ 
-k))C/ia(k, q), 




= (J, 






= A2(k)ty2Q(k, q) 






= (A2(k) - A3(q- 


-k))C/2a(k, q), 




= 0, 






= 0, 




T^f (k,q) 


= A3(k)C/|Jk,q) 





where we have used 
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C/2.(k,q)--[/,Vq-k,q) 



(4.12) 
(4.13) 



B. Mode coupling equations 

Using the standard mode couphng formahsm for the 
spin variables s" (q) one gets the fonowing equations for 
the diagonal elements of the memory matrix 



ra(q,i) = 



1 



(^^W,^^(o)) 



4J2 



T^(k,q)^^rfV(k',q) 

x(^kW^^^kW>4(0)<-k'(0))- (4.14) 

The off-diagonal elements are zero since the correspond- 
ing relaxation functions vanish due to the symmetry 
properties of the Hamiltonian. If the four-point corre- 
lation function on the right hand side is factorized into a 
product of two two-point correlation functions, one gets 



r„(q,i) = iMZ! / ^i^f^(k,q) 
Xa(q) 7k ^ 



/37 



x$^(k, i)$^(q-k,t), 



(4.15) 



with the vertex functions iff'' (k, q) for the decay of the 
mode a into the modes (3 and 7 given by 

Xf'^(k,q) = rr(k,q)C/f^(k,q) (A^(k)- A^(q-k)) , 
iff^(k,q)=rf''(k,q)rf^(k,q), /3 ^ 7, (4.16) 

The corresponding equation for the Fourier transform 
reads 

/oc 
dte''^^ r„(q, t)dt 
-00 



/37 



X 

where — J duj/2'K and 



(4.17) 



*a(q,t^) 



»Xa(q) 



(4.18) 



denotes the half-sided Fourier-transform of $q. (q, t) . If 
the transport coefficients vary only slowly with lo one may 
replace the relaxation functions by simple Lorentzians 



«'a(q,t^) = 



»Xa(q) 



(4.19) 



i.e., the transport coefficients are replaced by their values 
at cj = 0: 



ra(q) =r„(q,c^-0). 



(4.20) 



This additional approximation finally leads to a simpli- 
fied set of coupled integral equations for the diagonal 
elements of the transport coefficients at zero frequency: 



r,(q)=4fcBTj2^Xf^(k,q) 



X/3(k)x7(q-k) 



Xa(q) 



r/3(k) 



r^(q-k)- 



(4.21) 



As emphasized before, the uniaxial anisotropy and the 
dipolar interaction introduce two extra length scales 
and respectively. This entails the following exten- 
sion of the static scaling law for the spin susceptibilty (in 
the eigenvector basis) 



Xa(q;^,'7D,9A) = V^Xa(q;R-), 



1 

J? 



(4.22) 



where q is a unit vector and R — {x, y, z) a vector of 
scaling variables 



1 



y = 



<1l 



and 



(4.23) 



The suc cpti bil ites are given in terms of the eigenvalues 
in Eqs. (|3j)-(|lo|) as 



Xa(q;R) = 9^A„^(q;C,gD,g^). 



(4.24) 



Note that the susceptibilities of the spin variables Sq de- 
pend not only on the three scaling variables a;, y and z, 
but also on the direction q = q/g of the wave vector q 
with respect to the easy axis of magnetization. Thus, all 
scaling functions depend on Jour variables. The mode 
coupling equations are consistent with the same type of 
scaling law for the line widths rc((q; ^, (Zj^ , (7^) 



(4.25) 



where the dependence of the line width on the length 
scales has now been made explicit in its argument. In- 
serting the above scaling laws, Eq.( 4.22| ) and Eq.(4.25) 
and the scaling of the vertex functions 
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i^f^(k,q;C,<7,,gJ=g^i^^(p,q;R) 



(4.26) 



we get 



Jo Jo Jo Xa(q;-R) 



p5/2^^(r;,R/p) +p^_/S^(r;_,R/p_) 



(4.27) 



where p = k/q, p_ = (k — q)/g, and the first argument 
in the scaling functions for the hnc width now gives the 
azimuthal angle between the wave vector and the easy 
axis of magnetization. The non universal scale factor 



Ar2(27r)3 ' 



(4.28) 



and the dynamic exponent for the isotropic Heisenberg 
ferromagnet 

z = \. (4.29) 
Here we have introduced polar coordinates such that 

(4.30) 



cos /i sm V 
q = I sin ji sin v 

cos,v 



p cos <T sm rj 
p — \ p sin cr sin rj | , 

p cos J] 

/9_ COS (T- sm 77_ 
p = I /?_ sincT- sin?]- 
p_ cos r/_ 



(4.31) 

(4.32) 
(4.33) 



with 



P- = \/l + p^ ~ 2p{suiiysmr]cos{p. — a) + cos cos 77), 
sin p sin v — p sin a sin 77 



cr_ = arctan 



?7_ = arccos 



cos p sm V — p cos cr sm 77 
cos V — p cos 77 



C. Solution of the mode coupling equations 

Now we are going to discuss the solution of the mode 
coupling equations derived in the preceding section. 
Since the general case including uniaxial anisotropy as 
well as dipolar interaction is quite complicated we will 
first shortly review the limiting cases of a) uniaxial and 
b) dipolar anisotropy, respectively. Since both cases are 
discussed in detail in the literature we will be rather 
brief and mention only those aspects which will be rele- 
vant for the subsequent discussion. 



a. Uniaxial limit: The equations of motion in the 
uniaxial limit reduce to 



= -2 J{q' -2q-k~m) (4-34) 



-2 / J(g^-2q.k-77i)^£5^_k, 



5^ = -2 / J(g^-2q.k)5^5f 



q-k- 



(4.35) 
(4.36) 



The eigenmode analysis shows that the mode coupling 
equations (in Lorcntzian approximation) become diago- 
nal in the spin fluctuations Sa perpendicular and parallel 
to the easy axis of magnetization (a =_L, ||) 



r±(q) 
rii(q) 



AksTJ^ f (g2 - 2q • k - t 



-2\2 



XII (q) A rx(k)+r||(q-k) 
AkeTJ^ f (92-2q-k)2 



(q) Ar±(k) + ri(q-k)- 



(4.37) 
(4.38) 



The equations ( 4.37 ) and ( 4.38 ) obey the dynamical scal- 
ing laws [plpl 



(4.39) 



with the following mode coupling equations for the scal- 
ing functions of the line widths 



"/±{x,z) — 2tt / dp drj sin rj K±{p,ri;x,z) 
Jo Jo 

+ Pl 7||( — >— ) 
P P P- P- 



7|| (a;, z) = 27r / dp drj sin rj K\\{p,rj; x, z) 
Jo Jo 



,x z . 



P^ 7-l(-,-) + P- 7±( — , — ] 
P P P- P- 



and the scaled vertices, 

X^{p,x,z)x\\{p-,x) 



(4.40) 



(4.41) 



Ki_{p,'q]x,z) 



K\\{p,ri;x,z) 



X±{'^,x,z) 
Xp^ (1 - 2/9C0S77 - z^f, 

X±{p,x,z)x±{p-,x) 
X\\{'^,x,z) 

xp^ (1 - 2/9C0S77)^, 



(4.42) 



(4.43) 



the dynamic exponent z = 5/2 and the non- universal 
scale 



A 



lUkBTV^ 
(27r)3iV2 ■ 



(4.44) 



The crossover of the critical dynamic exponent is con- 
tained in the scaling functions 7a (2;, z). The above mode 
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coupling equations are essentially the same as those de- 
rived in Ref. |22| and their numerical results are con- 
firmed by our analysis here (see below). 

As summarized in table || the mode coupling equa- 
tions can be solved analytically in the uniaxial (U) and 
isotropic (I) critical (C) and hydrodynamic (H) limiting 
regions. These are defined by UC: z ^ 1, a; ^ 1; IC: 
z < 1, x < 1; UH: z > a;, a; > 1; IH: z < a;, a; > 1. 

To examine the uniaxial crossover precisely at the crit- 
ical temperature, Fig. ^ and Fig. ^ display the results of 
the mode coupling theory for the scaling functions for 
the spin fluctuations parallel and perpendicular to the 
easy axis of magnetization for T = Tc against the wave 
vector; i.e. z^^ — q/q^- We define an effective dynamic 
exponent by 



10 II I 

9 - 

8 - 

7 - 

6 - 

a - 

4 - 
3-1 

2 - \ 

1 

' 

01 23456789 10 

1/z = q/qA 

FIG. 4. Scaling function for tlie line width of tlie spin 
fluctuations perpendicular to the easy axis of magnetization 
'y±(x,z) as a function of g/g^ at the critical temperature 
x = (T = T,). 



^a{q\i,qA 



(X q 



a 



(4.45) 



The results of the mode coupling theory clearly show that 
the crossover from isotropic to uniaxial critical dynamics 
in the transverse line width (i.e. perpendicular to the easy 
axis of magnetization) occurs in the immediate vicinity of 
the uniaxial crossover wave vector q^ from the isotropic 
Heisenberg value Z(.g = 5/2 to Zcff — 0. In contrast, 
the crossover in the longitudinal line width (i.e. paral- 
lel to the easy axis of magnetization) from an effective 
exponent z^s — 5/2 to Zoff = 4 occurs at a wave num- 
ber larger than ^4 by approximately one order of magni- 
tude. It is interesting to compare this dynamical shift of 
the crossover position with the analogous problem of the 
crossover from isotropic Heisenberg to isotropic dipolar 
dynamics, where the shift of the crossover position in the 
critical mode (i.e. transverse to the wave vector of the 
spin fluctuations) is shifted to wave vectors smaller than 
the corresponding anisotropy scale, the dipolar wave vec- 
tor q^ . 




40 60 
1/z = q/qA 

FIG. 3. Scaling function for the line width of the spin fluc- 
tuations parallel to the easy axis of magnetization, 7|| {x, z), as 
a function of q/Qj^ at the critical temperature x — {T = Tc). 



The physical content of the two parameter scaling sur- 
faces is illustrated best by considering cuts for fixed q^ 
and various temperatures, since for a given material, 
q^ is fixed and the parametrisation by (p corresponds 
to a parametrisation in terms of the reduced tempera- 
ture {T - Tc)/Tc. In Figs. the scahng functions 
versus x — 1/q^ are displayed for different values of 
ip = arctan(g^^) = Ntt/20 with N = 0,1,.. .,9. For 
if = 0, corresponding to vanishing uniaxial anisotropy q^ , 
the scaling functions coincide with the Resibois-Piette 
scaling function If the strength of the uniaxial 
anisotropy is finite, the curves for the scaling func- 
tions approach the Resibois-Piette scaling function for 
small values of the scaling variable x and deviate there- 
from with increasing x. At fixed scaling variable x and 
with increasing temperature (increasing ip) the value of 
the scaling function for the spin fluctuations along the 
easy axis of magnetization is lowered with respect to 
the Resibois-Piette function, whereas it is increased for 
the scaling function of the spin fluctuations in the basal 
plane. 




FIG. 5. Scaling function for the linewidth of the spin fluc- 
tuations perpendicular to the easy axis of magnetization 7x 
as function of a; = 1/g^ for various angles ip = N 11/20, with 

= 0,...,9. 
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FIG. 6. Scaling function for the linewidth of the spin fluctu- 
ations parallel to the easy axis of magnetization 7|| as function 
of a; = 1/g^ for various angles ip — Nn/20, mit A'^ = 0, . . . , 9. 



b. Dipolar limit: As summarized in Table III the 
mode coupling equations in the pure dipolar limit can 
be solved analytically in the dipolar (D) and isotropic 
(I) critical (C) and hydrodynamic (H) limiting regions. 
These are defined by DC: y > 1, a: < 1; IC: y < 1, 
X < 1; DH: y > x, x > 1; IH: y < a;, a; > 1. 

Concerning the critical dynamical exponent one finds 
for the longitudinal line width a crossover from z = 5/2 
in the isotropic critical region to z = in the dipolar 
critical region, whereas for the transverse line width the 
crossover is from z = 5/2 to z = 2. The precise position 
of this crossover can only be determined numerically. 

A plot of the transverse and longitudinal scaling func- 
tions 7^(a:,2/) and ^^{x,y) can be found in Ref. For 
the dipolar crossover precisely at the Curie point the 
crossover from the isotropic Heisenberg to dipolar critical 
dynamics in the transverse line width occurs at a wave 
number, which is almost one order of magnitude smaller 
than the static crossover wave vector q^^ . The crossover 
of the longitudinal width, from z — 2.5 to z = 0, is more 
pronounced and occurs in the intermediate vicinity of 

. The reason for the different location of the dynamic 
crossover is mainly due to the fact that it is primarily the 
longitudinal static susceptibility which shows a crossover 
due to the dipolar interaction. Since the change in the 
static critical exponents is numerically small the trans- 
verse static susceptibility is nearly the same as for ferro- 
magnets without dipolar interaction. Hence the crossover 
in the transverse width is purely a dynamical crossover, 
whereas the crossover of the longitudinal width being 
proportional to the inverse static longitudinal suscepti- 
bility is enhanced by the static crossover. 

c. General case: For nonvanishing dipolar interaction 
and uniaxial anisotropy the mode coupling equations for 
the scaling functions of the spin fluctuations Sc((q) are 
given by 

/>00 pTT />27r 

7a(i';R) — dp drjsinri / da 
Jo Jo Jo 



where we have introduced polar coordinates 



(4.46) 



with 



R 



Re [0; 



R cos (fi sin 9 
R simp siii0 
R cos 9 



(4.47) 



^ G [0; 2 ] 



9 e [0; 



The mode coupling equations can be solved analytically 
only in some limiting cases. The results for these limit- 
ing cases are given table IV. In all other cases one has 
to solve the mode coupling equations numerically. For 
this it is best to write the scaling variables in spherical 
coordinates. The results for the three scaling functions 
are plotted in Figs. |^-||. In these figures we have plotted 
a table for each of the scaling functions, where in each 
element of the table the scaling function is shown as a 
function of the radial scaling variable R for a set of an- 
gles of the wave vector with respect to the z-axis. The 
three colums (from left to right) correspond to the angles 
(fi = tt/5, 27r/5 and 7t/2 and the rows (from top to bot- 
tom) correspond to 9 = Tr/5, 27r/5 and 7t/2. Note that 
ta.n ip = q^£, and tan 6* = ^^"^ +<ii/<lA- 

The continuous crossover from isotropic to uniaxial or 
dipolar behavior can then be easily read off from Figs. |^ 
I 

• For (p ^ Q and 9 ^ 71/2 all three scaling func- 
tions become identical to the Resibois-Piette scal- 
ing function for the isotropic Heisenberg ferromag- 
net. 

• In the uniaxial limit — > (corresponding to 
q^ 0) the first two scaling functions become 
identical. 

• Variation of the anisotropy strength q^ in the weak 
dipolar limit {q,j£, <C 1) (i.e., increasing 9 from 
to 7r/2 for small values of (p) gives the transition 
from the Resibois-Piette function to the uniaxial 
limit discussed above. 71 and 72 become identical 
to the scaling function for the uncritical transverse 
modes whereas 73 becomes the scaling function for 
the fluctuations along the easy axis. 



xp{p\ R) X7(p-;R) 
Xa(q; R) 



• For small uniaxial anisotropy (g^ <C ■^/^"■^ -I- or 
equivalently 9 close to 7r/2) one recovers the dipolar 
limit where 71 and 73 reduce to the two equivalent 
transverse modes, and 72 becomes the uncritical 
longitudinal mode. 

• Upon going along the diagonal from top-left to 
bottom-right in Figs. |^-^ one crosses over from the 
uniaxial to the dipolar limit and 71 and 73 become 
more and more similar. The scaling function 72 
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makes a continuous crossover from the scaling func- 
tion for the mode parallel to the easy axis (critical 
mode) in the uniaxial case to the scaling function 
of the longitudinal mode (uncritical mode) in the 
dipolar limit. 



• The first mode 71 does not show any dependence 
on the orientation v of the wave vector with respect 
to the z-axis. The reason is that this mode is by 
construction always perpendicular to both the z- 
axis and the wave vector q and hence idependent 
from the angle between these two axes. 



• The second and third mode show no i^-dependence 
neither in the dipolar nor in the uniaxial limit. The 
dependence on v is strongest close to 6* w 7r/4 and 
if « 7r/4. It is most pronounced for the third mode. 
This is due to the strong j/-dependenccof this mode 
in the hydrodynamic limit (see table IV). Whereas 
the scaling behavior of 71 and 72 does not depend 
on the orientation of the wave vector with respect 
to the 2;-axis in this limit. 

• The effective dynamic exponents Zoff become zero 
for all three modes in the case 7^ 0, except for 
73 in the case v = 7r/2. 



(p=7l/5 



(p=27l/5 



(p=7l/2 




0=71/5 



- 0=2ji/5 




0.0 I — ' — ' — ' — ' — ' — L 
0.0 2.0 4.0 6.0 

R 



FIG. 7. Scaling function 71 as a function of R for various angles 6 and tp. Each graph for fixed 9 and ip contains a set of 
curves parametrized by the angle v of the wave vector with respect to the z-axis: u — Nn/8 with A'^ = 0, ■ • • , 4. 
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(p=7l/5 



9=271/5 



(p=7i/2 



3.0 
2.0 
1.0 



0.0 



^^^^^^ 
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^ — fhi — ^ ^ 
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1,1 


^ h ' 
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^ — ) — ^ ^ 

1,1 



0=n/5 



e=2;i/5 



e=;t/2 



0.0 2.0 4.0 6.0 
R 



FIG. 8. Scaling function 72 as a function of R for various angles and cp. Each graph for fixed and contains a set of 
curves parametrized by the angle v of the wave vector with respect to the ^:-axis: v = Nn/8 with N = 0, - ■ ■ ,4. 



9=^/5 



(p=2;t/5 



9=71/2 




0=71/5 



6=271/5 



e=7t/2 



FIG. 9. Scaling function 73 as a function of R for various angles 6 and ip. Each graph for fixed 6 and ip contains a set of 
curves parametrized by the angle v of the wave vector with respect to the a-axis: v = Nn/8 with N = 0, - ■ ■ ,4. 
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d. Gd - no magnetocrystalline anisotropy: As we 
have already mentioned in the introduction Gd is a very 
interesting rare earth matrerial because it is an S-state 
ion with a large localized magnetic moment. Because 
of this it should have a very small magnetocrystalline 
anisotropy and be a much better model system for an 
isotropic Heisenberg ferromagnet than materials like Fe, 
EuO or EuS. Contrary to this expectation, however, ex- 
perimental observations teach us that Gd has an easy 
axis which coincides with the hexagonal axis of its hep 
structure. In the following section we will identify the 
dipole-dipole interaction in conjunction with the lattice 
structure as the source of this easy axis anisotropy. Then 
we can explicitely calculate the ratio of the two crossover 
wave vectors, 

qjq^= 7.8738, (4.48) 

and thus reduce the number of material parameters by 
one. This allows us to present our theoretical results for 
the line widths of the correlation functions in the same 
way as for the two limiting cases of the cubic dipolar and 
purely uniaxial systems. 

In Figures p^-p^ we compare the scaling functions for 
the line widths of the three modes 71 , 72 and 73 for the 
hep dipolar model with the corresponding results for the 
cubic dipolar model. Since the dipolar interaction is the 
dominant anisotropy the deviations of the results for the 
hep from the cubic model are small. For 71 and 73 we find 
that the relaxation times of the hep system are reduced 
as compared to the cubic model. This seems plausible 
since - as discussed in section III - the induced uniax- 
ial anisotropy supresses some of the critical fluctuations. 
This has to be contrasted with the behavior of 72, where 
the corresponding longitudinal mode in the cubic dipolar 
system is alread an uncritical mode. We find that the re- 
sult for 72 is almost identical to the longitudinal scaling 
function of the cubic dipolar system. 



T,.Yl 




0.5 

0.0 2.0 4.0 6.0 8.0 

x=1/q^ 



FIG. 11. Comparison of the scaling function 72 at = 
of the hep dipolar model (solid lines) with the longitudinal 
scaling function 7_l (dashed lines) of the cubic dipolar model. 
The scaling functions are given as a function of a; = l/g5 for 
various values of (fi — Ntt/20 with A'' = 0, 9. 




0.0 2.0 4.0 6.0 8.0 

x=1/q5 



FIG. 12. Comparison of the scaling function 73 at = of 
the hep dipolar model (solid lines) with the transverse scaling 
function jt (dashed lines) of the cubic dipolar model. The 
scaling functions are given as a function of a; = 1/q^ for vari- 
ous values of (p = Ntt/20 with N = 0, 9. 




0.5 

0.0 2.0 4.0 6.0 8.0 

x=1/q!; 



FIG. 10. Comparison of the scaling function 71 at = of 
the hep dipolar model (solid lines) with the transverse scaling 
function jt (dashed lines) of the cubic dipolar model. The 
scaling functions are given as a function of a:: = l/q(, for vari- 
ous values of (f = Ntt/20 with iV = 0, 9. 



The case v — tt /I is special since here the third mode 
in the hep model becomes critical (see section III). As a 
consequence there is no supression of critical fluctuations 
and one does not expect that dynamics is slowed down 
with respect to the cubic case. Indeed, as can be inferred 
from Fig. |l^, the mode becomes even faster than in the 
cubic case. 

The dependence on the relative orientation of the c- 
axis and the wave vector is not very pronounced. There 
is actually no i^-dependence for 71 since this mode is 
by construction always perpendicular to the c-axis. In 
Figs, p^p^ wc try to visualize the z^-dependence of 72 
and 73 at a set of parameters g^-,^ — 1.0 and = 1.0. 
These parameters do not correspond to Gd, because as 
we have already seen in Figs. the angle dependence 
is largest when and are comparable; hence it is 
quite small for the hep dipolar system Gd. 
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0.0 2.0 4.0 6.0 8.0 

x=1/q5 

FIG. 13. Comparison of the scaling function ^3 a.t v — ■k/2 
of the hep dipolar model (solid lines) with the transverse scal- 
ing function "fx (dashed lines) of the cubic dipolar model. The 
scaling functions are given as a function of a:: = for vari- 
ous values of If = Nn/20 with iV = 0, 9. 



Finally, let us discuss the behavior of the line width 
right at the critical temperature. Figs. p^-|l8| show the 
scaling functions of the three modes in the hep system 
compared with the corresponding results for the cubic 
dipolar and the uniaxial system. 

3.5 rr — ~ — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — I 



2.5 




0.5 I — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — I 
0.0 1.0 2.0 3.0 4.0 
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0.0 0.5 1.0 1.5 

x=1/q^ 

FIG. 14. Scaling function 72 for the hep dipolar model as a 
function of x — 1/q^ for fixed values qtj-,^ = 1.0 and — 1.0 
and a series of angles u — Nn/8 with A'' = 0, 1, • • • , 4 indicated 
in the graph. 
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0.0 1.0 2.0 3.0 4.0 

x=1/q5 

FIG. 15. Scaling function 73 for the hep dipolar model as a 
function of 2: = 1/g^ for fixed values = 1.0 and = 1.0 
and a series of angles v — Nn/8 with A'^ = 0, 1, • • ■ , 4 indicated 
in the graph. 



FIG. 16. Scaling function 71 at the critical temperature as 
a function oi q/q^^ (solid line) compared with the correspond- 
ing transverse scaling function -y^ of the cubic dipolar system 
(dot-dashed line) and the hard axis scaling function 7^ of the 
uniaxial system (dotted line). Note that in the uniaxial case 
the scaling function is plotted as a function of g/g^ but as a 
function of g/g^, for the dipolar cases. 




0.5 I ' ' ' 1 

0.0 1.0 2.0 3.0 4.0 

FIG. 17. Scaling function 72 at the critical temperature as 
a function of g/g^ for u — (solid line) and u — 71/2 (dashed 
line) compared with the corresponding longitudinal scaling 
function 7^ of the cubic dipolar system (dot-dashed line) and 
the hard axis scaling function 7^ of the uniaxial system (dot- 
ted line). Note that in the uniaxial case the scaling function 
is plotted as a function of g/g^ but as a function oi q/q^ for 
the dipolar cases. 



Again, we find that for the parameter values of Gd 
there is little difference between the results for the cubic 
and the hep dipolar system (but a huge difference to the 
results of a uniaxial model with no dipolar interaction). 
The most notable effect is that the crossover wave vector 
which marks the deviation from isotropic Heisenberg be- 
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havior (with a dynamic exponent z = 5/2) is shifted to 
larger values for the first mode and to smaller values for 
the second and third mode. Then there is also a clearly 
visible dependence on the orientation of the wave vector 
with respect to the c-axis for the third mode. 




FIG. 18. Scaling function 73 at the critical temperature 
as a function of q/qj^ for 1/ = (solid line) and v = n/2 
(dashed line) compared with the corresponding transverse 
scaling function 7^ of the cubic dipolar system (dot-dashed 
line) and the soft axis scaling function 7|| of the uniaxial sys- 
tem (dotted line). Note that in the uniaxial case the scaling 
function is plotted as a function of g/g^ but as a function of 
q/qj^ for the dipolar cases. 



V. COMPARISON WITH EXPERIMENTAL DATA 



The critical dynamics of Gd has been investigated ex- 
perimentally almost exclusively by several nuclear, i.e., 
hyperfine interaction (HFI), methods. The application of 
nuclear techniques to study critical phenomena in mag- 
nets has recently been reviewed by Hohenemser et al. . 
All of the hyperfine interaction methods are local probes 
which are related to a wave vector integral of the spin 
correlation functions. As such they offer a complement 
to neutron scattering. Dynamic studies, using hyperfine 
interaction probes, utilize the process of nuclear relax- 
ation produced by time-dependent hyperfine interaction 
fields, which refiect fluctuations of the surrounding elec- 
tronic magnetic moments. 

The hyperfine interaction stems from the magnetic in- 
teraction of the electrons with the magnetic field pro- 
duced by the nucleus. The hyperfine interaction of a 
nucleus with spin I, g-factor and mass with one 
of the surrounding electrons with spin S and orbital mo- 
mentum L can be written in the form Q 



hyp 



1 



I L 



Stt 



5(3) (x) I. S 



3(I-x)(S -x) 



(5.1) 



The first term represents the interaction of the orbital 
momentum of the electron with the nuclear magnetic 
moment of the nucleus. The second term is the Fermi 
contact interaction and the last two terms represent the 
dipolar interaction. The Fermi contact interaction is fi- 
nite only for electrons having a finite probability density 
at the nucleus, i.e. bound s-electrons or itinerant elec- 
trons. The Hamiltonian, Eq. (5.1), can also be used for 
the analysis of spin resonance experiments with muons 
(/uSR). However, these do not have bound electrons and 
hence the Fermi contact term involves only conduction 
electrons and is of the same order of magnitude as the 
(residual) dipolar interaction [p6[. 



A. /iSR measurement 

As shown in appendix ^ the muon damping rate can 
be written as 



A. 



q ^-q 



r<yi3r<vi 



$'^^(q), (5.2) 



where we have defined V — 7^(^o/47r)^(5L^B)^. The 
coupling of the muon spin and the spins of the magnet 
is described in terms of the coupling matrix G''^ , which 
reflects the particular symmetry of the lattice sites occu- 
pied by the muons. Since the most dominant contribu- 
tion to the damping rate comes from wave vectors close 
to the Brillouin zone center, the peculiar properties of the 
coupling tensor at small values of q will be impor- 
tant. The coupling tensor is determined by both the type 
of interaction between the muon spin and the spins of the 
magnet and the location of the muon in the lattice. As 
noted above the coupling contains dipolar interaction as 
well as a contribution from the Fermi contact field. Us- 
ing a decomposition into four orthorhombic sublattices 
{I = 0, 1,2,3) the dipolar interaction can be written in 
the form @, 



(5.3) 



where r^^; = i -I- r/ — Tq denotes the position on site i of 
sublattice I with respect to the position of the muon Tq. 
The lowest-order approximation in q of ^"^(q) is given 
by (see appendix ^ 



^"'^(q^ 0) = -Ait 



da 



(5.4) 



where we find for octahedral sites 



= dy = 0.3485 and 4 = 0.3030, 



and for tetrahedral sites 



d^ = dy = 0.3118 and 4 = 0.3764. 
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Adding the isotropic contribution from the Fermi contact 
field one finds in the limit q — > 



G 



Q/3 



-47r 
-At: 



qgqp 

q' 



with 



47r 



(5.5) 

Pa = rfa + _ " . (5.6) 

With the Fermi contact field Spc = -6.98 kG at T = K 
p6t , one gets n^H^/An = —0.278 |2^], and consequently 
for octahedral sites 

Px=Py = 0.0705 and = 0.0250, 

and for tetrahedral sites 

^ Py = 0.0338 and p^ = 0.0984. 

Before comparing the theoretical result with experiments 
one has to find a transformation between the reference 
frames used in the theoretical an alys is and the experi- 
mental setup, respectively. In Eq. 5.2 (i.e. in the experi- 
mental setup) one uses a reference frame (i, y, z) which is 



fixed with respect to the initial polarization of the muon 
beam, which is choosen to be along the z-axis. The refer- 
ence frame used in the theoretical analysis, {x, y, z) was 
chosen such that the z-axis coincides with the easy axis of 
magnetization, (see Fig. |l^). The transformation rules 
are given by 



X 

y 

z 



y cos a -1- z sma, 
—y sin a -I- z cos a. 



(5.7) 
(5.8) 
(5.9) 




z -« 




FIG. 19. Reference frame {x, y, z) for /iSR measurements, 
which is rotated by a around the a;-axis with respect to the 
reference frame {x,y,z). Here z is the direction of the easy 
axis of magnetization and z is the direction of the initial po- 
larization of the muon beam. 



Then the coupling tensor in the experimental reference frame reads 

,2 



G^^(q 
G**(q 
G*^(q 
G**(q 

G**(q 



0) 
0) 
0) 
0) 

0) 



-An 
-Att 
-Att 
-Att 

-An 



q2 



q-xqy , qxqz . 

— ^ cos a H 5— sm a— p^ 



■ sm a 



q- 

qxqy 



Qx qz 

— — cos a — Pa 



qxqz ■ 

COS a H — sm a 

q^ 



, COS a — Pz sm a 

„2 



iy 2 , c^qyq^ ■ , q'z ■ 2 

— COS a -I- 2-2^ sma COS a -I — ^ sm a — pycosa — pzSma 



G*^(q^O) = -47r 



ql ■ iy ■ , qyqz 

— sm a cos a 1 cos a sm a H — "hr- 

qi 



(cos^ a — sin^ a) — py cos a — Pz sin ( 



In the same way we find for the relation of the spin correlation functions in the reference frames defined 
to the initial muon beam polarization and the easy axis of magnetization, respectively, 

$=**(q) = ^'^''(q), 

$*?'(q) = $**(q) = $^2'(q)cosa + $='^(q)sina, 
(q) = (q) = (q) cos a ~ ^""v (q) sin a, 
$**(q) = $^^(q)cos^a-f $^^(q)sin^a + 2$^^(q)cosQ;sina, 

$6^(q-) = $^*(q) = <i>^''(q) sin a cos a - $'^^(q) sin a cos a + $^^(q)(cos^ a - sin^ a), 
(q) = (q) sin^ a -|- (q) cos^ a - 2$^^ (q) cos a sin a. 



(5.10) 
(5.11) 
(5.12) 
(5.13) 

(5.14) 

(5.15) 
with respect 

(5.16) 
(5.17) 
(5.18) 
(5.19) 
(5.20) 
(5.21) 



15 



In section [V C we have calculated the spin correlation 
functions ^""(q) in the eigenvector basis eQ,(q). Upon 
using the transformation, Eq. |4.6|, between cartesian spin 



components and the spin components in the eigenvector 
basis we find 



(5.22) 



With the Fluctuation-Dissipation-Thcorcm the spin cor- 
relation functions 'i>"'^(q) are given in terms of the stat- 
ical susceptibility and the linewidth, 



$""(q) 



'2kBT Xa(q) 
Mo(5lMb)' r"(q) 
2kBT 



JM9LtJ-B)^Aq-+^ A„(R) 7„(zy; R) ' 



(5.23) 



The muon relaxation rate Xz depends on the material 
parameters f o and characterizing the strength of 
the uniaxial anisotropy and the dipolar interaction, re- 
spectively. If we assume that the uniaxial anisotropy in 
Gd is solely due to the combined effect of dipolar inter- 
action and non cubic lattice structure, one can estimate 
the ratio of dipolar to uniaxial wave vector llq] 



— = 7.8738. 



(5.24) 



With this assumption the number of material parameters 
is reduced to one, g^^o- In comparing our theory with 
/iSR experiments at a polarization a — 90° we get the 
best fit to the data with g^^o = 0.13. This results in the 
following values for the uniaxial and dipolar wave vector 



q^ = 0.0165/Co, Qo = 0.13/^0. 



(5.25) 



The corresponding crossover temperatures, qa.dC — 1; 
are given by 



Ta = Tc + 0.43 K, Td=Tc + 16.54 K. 



(5.26) 



These set of parameters suggest the following crossover 
scenario. For T ^ To we expect critical behavior dom- 
inated by the (isotropic) Heisenberg fixed point. The 
relaxation rate shows power law behavior 



X (xt 



(5.27) 



with an exponent wi w i^(z — 1) w 1. For temperatures 
in the interval To > T > Ta dipolar interaction be- 
comes important. B ut, fro m our analysis of the uniaxial 
crossover in section IV C we have seen that the uniax- 



ial crossovers in dynamics sets in at wave vectors much 
larger than expected from an analysis of the static quanti- 
ties; i.e., the dynamic crossover in the longitudinal scal- 
ing function right at Tc is located at gcross ~ 10 x g^. 
Therefore, even for T > Ta we expect to observe effects 
from dipolar interaction as well as uniaxial anisotropy. 
Finally, for T < Ta the critical dynamics is determined 



by the uniaxial dipolar fixed point. Then the static sus- 
ceptibilities do no longer diverge for q ^ and T ^ Tc 
except when the wave vector q is perpendicular to the 
easy axis of magnetization. Since the relaxation rate 
is given by an integral over the whole Brillouin zone, 
the relative weight of the critical axis along which the 
susceptibility diverges becomes vanishingly small. As a 
consequence the relaxation rate Xz no longer diverges for 
T^Tc. 

Since the interaction between the spin of the muon 
and the lattice spins is a combination of Fermi contact 
and dipolar interaction, the muon relaxation rate is a 
rather complicated function of the relaxation rates along 
the eigendirections. Therefore, the temperature depen- 
dence of the relaxation rate can no longer be described 
in terms of simple power laws but shows a more compli- 
cated functional dependence. In Figs. ^ and |l] we show 
a comparison between the theoretical and experimental 
results for two different initial polarizations. 




T-Tc [K] 

FIG. 20. Experimental and theoretical results of the relax- 
ation rate A for tetrahedral and octahedral muon sites with 
a = 90°. 



octrahedral sites 
- tetrahedral sites 
measurement 89 - 




T-Tc [K] 

FIG. 21. Experimental and theoretical results of the relax- 
ation rate A for tetrahedral and octahedral muon sites with 
a = 0°. 



In Fig. 20 the initial polarization is inclined by an an- 
gle a — 90° with respect to the easy axis of magnetiza- 
tion. The solid line and dashed line are the theoretical 
result for the muon relaxation rate if the muons penetrat- 
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ing the sample are located at tetrahedral and octahedral 
interstitial sites, respectively. The comparison between 
theory and experiment favors tetrahedral sites. This is 
confirmed by /iSR experiments with the initial polariza- 
tion along the easy axis of magnetization (see Fig. pT] ). 
The ratio \^{90°)/ X^{0°) for T ^ becomes 1.2 and 
0.7 for octahedral and tetrahedral sites, respectively. The 
experiment is closer to the latter value. 



B. PAC and MS measurements 

Here we consider those hyperfine interaction probes 
(ME, PAC, NMR), where the Fermi contact term gives 
the dominant contribution to the hyperfine field at the 
nucleus. Then the corresponding interaction Hamilto- 
nian reduces to 



i7(i) = AFcrmil-S(t) 



(5.28) 



where I denotes the nuclear and S the electronic spin. In 
hyperfine interaction experiments one observes the nu- 
clear relaxation rate due to the surrounding fluctuating 
electronic magnetic moments. The standard experiments 
are performed in the motional narrowing regime, where 
the nuclear relaxation rate is directly proportional 
to the (averaged) spin autocorrelation time Tc 



-, poo -, 

2 J ^rf^3E'^""(^ = 0'*)' 



(5.29) 



where the spin autocorrelation function is given by 



$«"(R,t)-i({5"(x,t),^"(0,0)} 



(5.30) 



Hence the above hyperfine interaction methods provide 
an integral property of the spin-spin correlation hmction. 
Upon using the fluctuation dissipation theorem (FDT) 
we get 



Va 



Bz 3 4:^ ra(q;f,go,g^)' 



(5.31) 



The q-integration extends over the Brillouin zone (BZ), 
the volume of which is Vq. 

Important information about the behavior of the auto- 
correlation time can be gained from a scaling analy- 
sis. Upo n using the static and dynamic scaling laws 
Eq. (5.31) can be written as 



Tc (X Att / dqq 



Asr^ Xa (q; qt, q/ gp , g/g^ ) 
3 4:^ 7a {(i;q^,q/qo^q/qA) 



(5.32) 



where we have neglected the Fisher exponent rj. If there 
were no dipolar interaction and no uniaxial anisotropy, 
one could extrac t the temperature dependence from the 
integral in Eq. ( 5.32| ) with the result Tc oc This 
expression can be used to define an effective dynamical 



exponent Zcs(t), which depends on the correlation length 

by 



Tc OC CX 



T-Tc 
Tc 



(5.33) 



with Wcs = I'izcs - 1). 

If dipolar interaction and uniaxial anisotropy are ab- 
sent, one would expect to observe a critical exponent for 
the relaxation rate w « 0.70 x (5/2— 1) « 1.0. Dipolar in- 
teraction is known to be a relevant perturbation with re- 
spect to the Heisenberg fixed point; it leads to asymptotic 
static critical exponents which are only slightly differ- 
ent from the corresponding Heisenberg values, but since 
dipolar interaction implies a non-conserved order param- 
eter the asymptotic dynamic exponent becomes zd ~ 2. 
Hence dipolar interaction would induce a crossover from 
wi w 1.0 to wu ~ 0.7. Uniaxial interaction is also known 
to be a relevant perturbation with respect to the Heisen- 
berg fixed point. Again, the static critical exponents are 
not changed very much, e.g. one finds vi — 0.63, but the 
dynamic exponent becomes zj « 4 if the order parameter 
is conserved (zi « 2 otherwise). The corresponding ex- 
ponent for the hyperfine relaxation rate would turn out 
to be wi « 1.89 and wi « 0.63 for conserved and non- 
conserved order parameter, respectively. According to 
these scaling arguments it is hard to think of any dy- 
namic universality class which could lead to an effective 
exponent Wqs smaller than about 0.6. Actually, how- 
ever Mossbauer studies and PAC measurements on Gd 
show distinctly anomalous low values w ~ 0.5, which 
can not be explained by either of the above scenarios. 
This experimental puzzle can be resolved if one consid- 
ers the combined effect of dipolar interaction and uniax- 
ial anisotropy. As we have seen in our analysis of the 
static critical behavior of uniaxial dipolar ferromagnets, 
all the eigenvalues of the susceptibility matrix remain fi- 
nite upon approaching the critical temperature except 
when the wave vector of the spin fluctuations is perpen- 
dicular to the easy axis of magnetization. Since this is 
only a region of measure zero in the Brilloui n zon e one ac- 
tually expects that the relaxation rate Eq. ( 5.31 ) does no 
longer diverge upon approaching Tc, i.e., wud = 0. For a 
quantitative comparison with the experiment |^-|30|] we 
rewrite the auto correlation time in scaling form 



JVA 



1 



X I dv sin v — 

/o 3^ A„(i^;R)7"(i^;R) 



(5.34) 



where we have introduced polar coordinates. The lower 
cutoff is given by 



i?n = 



1 



S. + ql+ql, 



(5.35) 



where qsz is the boundary of the Brillouin zone. In 
the critical region it can be disregarded and replaced 
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by i?o = 0, since q^^ ^ 1a ^^"^ the integrand in 
Eq. ( 5.34| ) is proportional to ^/R for small R. For very 
small ^ (outside the critical region) the cutoff reduces the 
autocorrelation time with respect to the critical value. 
One should note, that the dominant wave vectors con- 



tributing to the relaxation time Tc in Eq. (5.34) are close 
to the zone center ||2^ . 

Let us now compare with hyperfine experiments on Gd 
mentioned above |p9|,^ ,^. The autocorrelation time Tc 
is shown in Figs. |22| and |23| for PAC experiments and 
Mossbauer spectroscopy, respectively. Both set of data 
are in good agreement with the results from mode cou- 
pling theory for T - Tc < lOK. Note that besides the 
overall frequency scale there is no fit-parameter, since we 
have used the same set of values for the dipolar and uni- 
axial wave vector as for our comparison with /iSR exper- 
iments. At higher temperatures the PAC data are above 
and the Mossbauer data below the theoretical prediction. 























Td 


\ 



0.0001 0.01 0.1 1 10 100 1000 

T-Tj- [K] 

FIG. 22. Experimental and theoretical results of the auto- 
correlation time Tc for PAC experiments. 




T-Tc[K] 

FIG. 23. Experimental and theoretical results of the auto- 
correlation time Tc for Mossbauer spectroscopy experiments. 

The dipolar and uniaxial crossover temperatures indi- 
cated in the figures show that all the data are within the 
crossover regime between isotropic dipolar and uniaxial 
dipolar critical behavior. This explains the anomalous 
low value for the effective exponent w ~ 0.5. This re- 
sult clearly shows for the first time that the universality 
class for the critical behavior of Gd is the uniaxial dipolar 



ferromagnet. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the critical dynamics of three- 
dimensional ferromagnets with uniaxial anisotropy tak- 
ing into account exchange and dipole-dipole interaction. 
This analysis was mainly motivated by the puzzling ex- 
perimental situation for the rare earth material Gd. Due 
to the fact that it is an S-state ion one would expect that 
it should have a vanishing magnetocrystalline anisotropy 
and therefore be an almost ideal model system for an 
isotropic Heisenberg ferromagnet (model J). But, the 
actual experimental observation shows large deviations 
from the expected behavior. As explained in detail in 
the main text, the way out of this puzzle is to realize 
that there is always a dipolar interaction between the 
magnetic moments in the material which in conjunction 
with the crystalline anisotropy can lead to anisotropy in 
the magnetic properties of the system. 

The dipolar interaction and uniaxial anisotropy in- 
troduce two length scales, the dipolar wave vector 
and the wave vector measuring the uniaxial 
anisotropy. Note that there are two sources for the uniax- 
ial anisotropy: magnetocrystalline (non S-state character 
of the magnetic ions) and dipolar. For Gd there is only 
the latter and hence q^ and are not really indepen- 
dent but their ratio is fixed, q^^/qj^ = 7.87. Both length 
scales then have a common source, the dipole-dipole in- 
teraction. 

The presence of these length scales leads to a quite 
complicated crossover scenario already for the static spin- 
spin correlation function. The essence of the physics, 
however, can already be seen within an Ornstein-Zernike 
approximation. Without uniaxial anisotropy two out of 
the three eigenmodes of the correlation function are criti- 
cal, whereas the third correlation function starts to satu- 
rate once the wave vector becomes smaller than the dipo- 
lar wave vector . It is quite remarkable that due to 
the combined effect of dipolar and uniaxial anisotropy 
all eigenvalues of the susceptibilty matrix remain finite 
in the long wave length limit and upon approaching the 
critical temperature. Only if the angle v between the easy 
axis of magnetization and the wave vector is v — 90° the 
third eigenvalue (corresponding to the transverse mode 
when uniaxiality is neglected) becomes critical. 

We have used mode coupling theory to derive a set of 
integral equations for the time-dependent spin-spin corre- 
lation functions. Using the Lorentzian approximation we 
have solved these equations numerically and determined 
the complete dynamic crossover scenario. Spezializing to 
Gd with a relatively weak uniaxial anisotropy (induced 
by the dipolar interaction) we find that the scaling func- 
tions for the line widths show a behavior which is similar 
to what is found for cubic dipolar systems (with no uni- 
axial anisotropy). The deviations become largest as the 
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critical temperature is approached. 

In the final section we have compared our results with 
/iSR, PAC and Mossbauer experiments and find quan- 
titative agreement with our theory. This explains the 
anomously low values of the critical exponents found in 
these experiments as a combined efi^ect of lattice structure 
and dipolar interaction. From the quantitative agree- 
ment between theory and experiment the following con- 
clusions can be drawn: 

(i) The universality class of Gd is the uniaxial dipolar 
ferromagnet. 

(ii) The dominant factor for the uniaxial anisotropy in 
Gd is solely the dipolar interaction. 

(iii) The theory even allows us to predict that muons 
in Gd are located at tetrahedral interstitial sites close to 
Tc- We expect that the analysis presented in this paper 
will be valuable for interpreting dynamic measurements 
on the non-cubic magnetic system. 



with the dipolar contribution 



3r? 



(A3) 



and the Fermi contact contribution H^^ , which equals 
H for nearest-neighbor atoms and zero otherwise. The 
number of nearest neighbors equals n^. ^'^^ is the sym- 
metrized spin-spin correlation function between spins at 
sites ii and 12 and zero frequency. 



<l = = Q)Sl) + {Sl{.^ - 0)5f^(c^ = 0))] (A4) 

Now we are going to rewrite Eq. ( |Al[ ) in Fourier space. 
This leads to the evaluation of dipole sums 



1, 



(A5) 
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APPENDIX A: ZERO-FIELD ^iSR 
DEPOLARIZATION RATE AND SPIN 
CORRELATION FUNCTIONS 



How one deals with such dipole sums for cubic lattices 
has been described in Ref. . Here we have to calculate 
the dipole sum for lattices with a hep structure. This is 
most conveniently done by decomposing the lattice into 
four orthorombic sublattices such that each site is in ex- 
actly one sublattice ||l^ . The basis lattice vectors of the 
sublattices in terms of the basis vectors of the hep lattice 
are given by 



ai = a, a2 = 2b + a, and aa — c. 



(A6) 



The volume of an orthorhombic unit cell fortho is four 
times the volume Vf^cp per atom in gadolinium, 



1 



rWortho- 



(A7) 



It can be shown [^l| that the muon damping rate in 
longitudinal geometry is related to the spin-spin correla- 
tion function of the host material by a sum over all lattice 
sites 



y2 



G-^(r,JG-'^(r,.Jci>f; 



13-1 
i2 



ili-2 13-t 



/37 

42 



(Al) 



where is the vector pointing from the interstitial po- 
sition of the muon to the lattice site i. Here we have 
defined V = {^iol^TrY"ff^{gLlJ'BY ■ Note that the sum is 
weighted by the coupling tensor G"'^ (r^ ) of the muon and 
the spins of the magnet, 



G"^(r,) = i?"''(r,)+n^i?(r.)5o/3, 
Putting things togehter we find 



(A2) 



A convenient choice of the origins of the four sublattices 
in terms of the basis vectors ai , a2 and aa is 



Rn = 



Ri = -ai 
2 



-a2 

2 



R2 
R3 



—sl-2 H — a^ 
3 2-^ 



7^1 



1 



:a3. 



2 6 2 
The muon can either be located at a tetrahedral site 



(AS) 
(A9) 

(AlO) 

(All) 



or a octahedral site 



To = ^ai + -a.2 + Tag. 
2 D 4 



(A12) 



(A13) 
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= ^EEE [G^''(ro + ii + R;JG-'^(ro + i2+R;.) + G^''(ro + ii + R;JG^^(ro + i2 + R;J] (1-14) 
Introducing Fourier transformed spin variables for each of the sublattices 

5"(ro + i + RO = ^ E ^^P[^q -(^0+1 + R;)] ^z"(q) 
q 

where R; are the lattice vectors of the Z-th sublattice, one gets 
2 ^^^o\^ {qlI^bY 



(1.15) 



[G^''(q)G^^(-q) + G'^'^(q)G'^^(-q)] $'='^(q) 

q p-f 



(1.16) 



where the Fourier transform for the symmetrized spin-spin correlation function and the coupling tensor reads 

3 



G"^(q) = E E ^"'"(^o + i + RO exp[iq • (ro + i + R;)] 



1=0 i 



/37 



^ ^ exp[zq • (ii - ia)] exp[iq • (R,, - R, J]$'3^(q) 
q 

^ J2 exp[^q • ((ro + R/^ + ii) - (ro + R;, + i2))]$''^(q) 



iV2 

q 

^''"(q) = ^ [{S^i^ = 0)5Zq(c. = 0)) + {Sl^iu; = 0)S^{co = 0))] 



(1.17) 
(1.18) 
(1.19) 
(1.20) 



The summation over q runs over the N vectors in the in /iSR measurements close to Tc the dominant contribu- 
first Brillouin zone. For sufficiently large N the summa- tion to A comes from the vicinity of the Brillouin zone 

center. The coupling tensor G°'^{q) can therefore be ap- 
proximated by its limiting behavior near q 0. 



tion may be replaced by an integral 

= ^ / E K'G^-\ + Gf G^JJ <I>^^(q). (1.21) 



Now we turn to the evaluation of the dipolar sums in 
the coupling tensor. The dipolar sums can be evaluated 
In principle the damping rate depends on the wave in Fourier space by the method of Ewald summation, 
vector dependence of the coupling tensor and the sym- The sum in Fourier space is divided in a sum over the 
metrized spin-spin correlation function over the whole direct lattice and a sum over the indirect lattice, so that 
Brillouin zone. It was shown, however, in Ref. that both sums converge quickly. 

With the decomposition into four orthorombic sublattices we get 

3 



^"''(q) = E E exp [zq • (r,,,)] 

;=o i 

3 

= 51 E ^""(^ + ^o) • (i + R; - i-o)] 



1=0 i 

3 

Vhcp ^ exp[iq • (ro + R/)] 

;=o 



9^ / exp(iq • i) 
dxa dxp \^ |i - x| 



/ J x=ro-Ri 

The latter expression can be evaluated in the same way as for cubic lattices p3] . One finds |31[ 



1=0 



1=0 



= -47r 



-G"'^(q) 



(1.22) 



(1.23) 



where 
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L^^K^ + q^)^Kp + qp)ipo (^^^^) exp[-iK • (R^ - ro)] 



4p2 



3 

[^P'i^l + i - ro)a(R; + i - ro)/3 <^3/2(p^r-2) _ (^./^(pV^) exp[iq • (R; + i) - ro]] . (1.24) 



To lowest order in q this reduces to 




9a 9/3 



- C"^(q = 0) 



0.3485 
0.3485 
0.3030 



0.3118 

0.3118 
0.3764 



(1.25) 



(1.26) 



(1.27) 
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Pi — P3 



0.44 



0.53 



-0.16 
0.44 
4.12 
-0.06 



-0.11 
0.31 
4.32 



TABLE I. Coefficients in the expansion of the dipolar ten- 
sor for a hep lattice structure with c = 5.78A, a — 3.62^4, 
and c/a = 1.59, as app ropriate for Gd bear Tc- Values are 
obtained from Ref. Ilq] upon multipying with the volume of 
the primitive cell Va- 



region 








DC 
DH 
IC 


,5/2 




,5/2 


IH 


^-1/2,2 




^-1/2,2 




C : 


critical 


D : 


dipolar 


H : 


hydrodynamic 


1 : 


isotropic 



TABLE III. Analytic results for the asymptotic behavior of 
the transverse and longitudinal linewidths for isotropic dipo- 
lar ferromagnets, where the spins are coupled by short-range 
exchange interaction and long-range dipolar interaction. 



region Fx T 



uc 








UH 








IC 


,5/2 




,5/2 


IH 


^-1/2,2 




^-1/2,2 




C : 


critical 


U : 


uniaxial 


H : 


hydrodynamic 


1 : 


isotropic 



TABLE 11. Analytic results for the asymptotic behavior 
of the transverse and longitudinal line widths for aniotropic 
ferromagnets with a uniaxial magneto-crystalline anisotropy, 
where the spins are coupled by short range exchange interac- 
tion only. 



region 


Li 


r2 


Ts 


UDC V / 90° 




iT 




UDC V = 90° 




iT 


1T1' 


DUG V / 90° 




iT 


iT 


DUG V = 90° 




iT 




UDH V / 90° 




iT 


iT 


UDH u = 90° 




iT 




DUH V / 90° 




iT 




DUH V = 90° 




iT 




UC 


iT 


iT 


i-/^'i' 


UH 


iT 


iT 




DC 




iT 




DH 


1/2 f-2 


iT 




IC 


,5/2 


,5/2 


,5/2 


IH 


^-1/2,2 


^-1/2,2 


^-1/2,2 



C : critical D : dipolar H : hydrodynamic 
U : uniaxial I : isotropic 

TABLE IV. Analytic results for the asymptotic behavior 
of the linewidths in the eigendirections for ferromagnets with 
a uniaxial magneto-crystalline anisotropy, where the spins are 
coupled by short-range exchange interaction and long-range 
dipolar interaction. 
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